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a^ : 

5^ ■ ABSTRACT 

H : 

u ■ The BL Lacertae object Markarian 501 was identified as a source of 7-ray emission 

^ ■ at the Whipple Observatory in March 1995. Here we present a flux variability 

. analysis on several times-scales of the 233 hour data set accumulated over 213 nights 

, (from March 1995 to July 1998) with the Whipple Observatory 10 m atmospheric 

I Cerenkov imaging telescope. In 1995, with the exception of a single night, the flux 

QQ I from Markarian 501 was constant on daily and monthly time-scales and had an 

' average flux of only 10% that of the Crab Nebula, making it the weakest VHE source 

• detected to date. In 1996, the average flux was approximately twice the 1995 flux and 

0^ ■ showed significant month-to-month variability. No significant day-scale variations were 

. detected. The average 7-ray flux above ~350 GeV in the 1997 observing season rose 

Q^. to 1.4 times that of the Crab Nebula - 14 times the 1995 discovery level - allowing 

Q I a search for variability on time-scales shorter than one day. Significant hour-scale 

I variability was present in the 1997 data, with the shortest, observed on MJD 50607, 

^ I having a doubling time of ~2 hours. In 1998 the average emission level decreased 

^ ' considerably from that of 1997 (to ~20% of the Crab Nebula flux) but two significant 

■ flaring events were observed. Thus, the emission from Markarian 501 shows large 

. amplitude and rapid flux variability at very high energies as does Markarian 421. It 

. , also shows large mean flux level variations on year-to-year time-scales, behaviour which 
has not been seen from Markarian 421 so far. 
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Subject headings: BL Lacertae objects: individual (Markarian 501) — gamma-rays: 
observations 



Introduction 



Three active galactic nuclei (AGN) have been discovered to be very high energy (VHE, E^300 
GeV) 7-ray sources by the Whipple Observatory 7-ray collaboration: Markarian 421 (Mrk 421) 
(punch et al. 1992|) , Markarian 501 (Mrk 501) (|Quinn et al. 1996|) and lES 2344+514 (|Catanese et| 
al. 199§| ). These are the three closest BL Lacertae objects (BL Lacs), with redshifts in the range 
0.0308 - 0.044, and are among the brightest at X-ray energies. They are all classified as X-ray 
selected BL Lacs (XBLs) as their synchrotron spectra extend into the X-ray range. A fourth 
BL Lac, PKS2155-304, has been detected in VHE 7-rays by the University of Durham group 
(IChadwick et al. 1999| ). 



The Energetic Gamma-Ray Experiment Telescope (EGRET) on board the Compton 
Gamma-Ray Observatory has detected at least 51 AGN at energies >100 MeV (Thompson et al, 



1995t pVIukherjee et al. 1997] ). They are all members of the blazar class of AGN, which include 



flat spectrum radio quasars and BL Lacs. Of the EGRET-detected blazars, 14 are BL Lacs with 
12 Radio Selected BL Lacs (RBLs) and only two XBLs. Mrk 421 is the only VHE source in 
this catalogue and it is among the weakest. However, Mrk 501 has recently been detected at 
the 4cr level with EGRET ( [Kataoka et al. 1999 ). One of the most striking characteristics of the 



EGRET-detected blazars is variability; 42 of the 51 AGN exhibit variability ( Mukherjee et al 



1997). Variability time-scales as short as 4 hours have been observed (Mattox et al. 1997) 



VHE observations of Mrk 421 and Mrk 501 have also revealed extreme variability (e.g. Gaidos 



et al. 1996; Quinn et al. 1996). The VHE flux has been measured to vary by nearly a factor of 



100 in Mrk 421 (|McEnery et al. 1999| ) and, as we show in the following sections, Mrk 501 has 
been measured with fluxes ranging from 0.1 to 5 times the Grab Nebula flux with the Whipple 
Observatory atmospheric Cerenkov telescope. The large collection area (~ 3.5 x 10^ m^) of the 
Whipple Observatory telescope permits sensitive studies of variability on time-scales inaccessible 
to space-based telescopes. Indeed, the shortest observed variability of any blazar at any 7-ray 



energy, a 30 minute duration flare observed from Mrk 421 ( Gaidos et al. 1996 ), was measured 
with this telescope. Both Mrk 421 and Mrk 501 have been closely monitored by the Whipple 
GoUaboration since their discovery, with an ~0.5 hr exposure/night being sufficient for detection 
of flaring activity. Prior to 1997, the 7-ray emission from Mrk 421 was generally observed to have 
a higher mean flux ( pchubnell et al. 1996 ) and to have been more frequently variable (Buckley et 
al. 19961 ) than that of Mrk 501. 



In Spring 1997, the Whipple GoUaboration observed Mrk 501 to be in an unprecedented 
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high emission state at VHE energies, as subsequently confirmed by several independent Cerenkov 
imaging groups ( Protheroe et al. 199^ ). A public Compton Gamma-Ray Observatory Target 



of Opportunity was initiated in response to a request by the Whipple Collaboration in 1997. 
Evidence of correlated variability in data from the Whipple Observatory 7-ray telescope, the 
Oriented Scintillation Spectrometer Experiment (OSSE) and the All-Sky Monitor (ASM) of 
the Rossi X-ray Timing Explorer is presented in Catanese et al. (1997). During this observing 
campaign the energy output of Mrk 501 in VHE 7-rays was comparable to that in the 2-100 keV 
range but the variability amplitude was larger. The correlations seen may imply some relativistic 



beaming of the emission, given that the spectrum extends to ;^7 TeV ( Samuelson et al. 1998 ). 
There was also some indication that the optical U-band flux was higher on average in the month 
of peak 7-ray activity. 

Here, based on 4 years of data, we present a study of the variability of the 7-ray flux above 
~350 GeV from Mrk 501. Details of the observations are given in ^ and the analysis methodology, 
including 7-ray selection criteria and the test for variability, is described in §^. The results of the 
analysis are presented in §Q and their implications briefly discussed in 



2. Source Observations 

Observations were made with the Whipple Observatory 10 m atmospheric Cerenkov imaging 
telescope ( Pawley et al. 1991 ), located on Mt. Hopkins in southern Arizona. The 10 m reflector 



images the Cerenkov radiation from cosmic-ray and 7-ray initiated air-showers onto a high 
resolution camera mounted in the focal plane. Subsequent off-line analysis of the images (described 
below) facilitates the selection of candidate 7-ray events. 

The high resolution camera utilises fast photomultiplier tubes (PMTs) arranged in a 
hexagonal array, with inter-tube spacing of 0?25. During 1995 and 1996 the camera consisted 
of 109 PMTs with a resulting field of view (FOV) of ~ 2? 8. An event was recorded when any 
two of the inner 91 PMTs registered a signal ;^40 photoelectrons within an effective resolving 
time of 15 ns. For the 1997 observations a further 42 PMTs were added, resulting in a FOV 
of ~ 3?4. The trigger condition remained the same as for the 109 PMT camera. In Summer 
1997 a new camera, containing 331 pixels, was installed. This camera has a FOV of ~ 4?8. The 
trigger condition for this enlarged camera required that any two of the 331 pixels produce a signal 
;^40 photo-electrons within an effective resolving time of 8 ns. The telescope was also triggered 
artificially, after every 24 events for the 109/151 pixel cameras and once every second for the 331 
pixel camera, to determine the background sky-brightness level in each PMT. Light-cones, which 
minimise the dead-space between PMTs and reduce the albedo effect, were used on the 109 and 
151 pixel cameras, but were not yet installed on the 331 pixel camera. 

In general, two modes of observation are used: on/off and tracking. With the on/off 
mode the source is tracked continuously for 28 minutes and then, to estimate the background, a 
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region offset in right ascension (RA) by 30 minutes (allowing 2 minutes slew time) is tracked. This 
has the disadvantage that an equivalent amount of observation time is spent looking away from 
the source. Alternatively, the tracking mode, where the background is estimated from the ON 
source run itself (see § |3.2D , can be used. In this case an off source run is not required for each 
ON source run, allowing continuous monitoring of an object. However, off runs are still needed to 
determine the response of the telescope to background events. 

Observations are typically made when the source zenith angle is less than 35° and are referred 
to as small zenith angle (SZA) observations. Observations at large zenith angles (LZA, typically 
55° to 70°) may also be made. Increasing the zenith angle has the effect of increasing the energy 
threshold, but has the benefit of increasing the collection area. Thus, it is an excellent method 
for increasing photon statistics to facilitate the determination of the energy spectrum at higher 
energies. For a detailed description of the LZA technique see Krennrich et al. (1997). 



Since its discovery as a 7-ray source in 1995 (g uinn et al. 1996|), the VHE 7-ray emission 



from Mrk 501 has been monitored intensively with the Whipple Observatory 10 m telescope. 
Only SZA observations taken under good sky conditions have been considered in the analysis for 
variability presented here. Our selection includes data from 56, 50, 55 and 49 nights of observation 
in the Spring - Summer periods of 1995, 1996, 1997 and 1998 respectively. Table |l] summarises 
the resulting database. 



3. Data Analysis 

3.1. 7-ray Selection 

The vast majority of events detected by Cerenkov telescopes are cosmic rays. Candidate 
7-ray events are selected on the basis of the shape and orientation of the Cerenkov images: 7-ray 
images are typically more compact and elliptical than background hadronic images and tend to 
have their major axes aligned with the source location in the FOV. Background cosmic ray events, 
on the other hand, have random orientations. 



Each image is first subjected to a cleaning procedure ( Fegan 1997] ) which suppresses pixels 



which are dominated by light from fluctuations of the night-sky background. A moment-fitting 
routine is then used to calculate various image parameters. The shape of each image is 
characterised by the parameters length and width and the orientation by a, the angle between 
the major axis of the image and the line joining the source location in the FOV to the centroid 
of the image. In addition, for the data taken in 1998, an asymmetry cut has been included. 7-ray 
images have a cometary shape with a tail which points away from the source location in the FOV 



( Buckley et al. 199S ) and thus their intensity profiles have positive asymmetry. The larger FOV 
of the camera used in 1998 allows this parameter to be accurately determined, something which 
was not possible with the smaller FOV cameras. 



- 5 - 



Before the application of shape and orientation cuts a software trigger cut is also applied to 
eliminate events close to threshold, some of which are induced by noise fluctuations. The software 
trigger involves cuts on the image size (i.e. the total number of photoelectrons recorded), the 
counts in each of the brightest two tubes {maxl, max2)^ and a requirement that at least three 
tubes above a low noise threshold (2.25(T, where a is the RMS sky-noise in a PMT, as determined 
from the artificially triggered events) be neighbours {NBR3). A distance cut is applied to eliminate 
images which are too close to the center of the camera and will have poor a reconstruction and 
also those events which have occurred too close to the edge of the FOV and may be truncated. 

Due to the continuous evolution of the high resolution camera, and changes such as 
deterioration of mirror reflectivity through weathering and the presence or absence of light-cones, 
the optimum data analysis cuts differ for each year. The cuts used for a given telescope 
configuration are optimised on an independent data set, usually data taken on the Crab Nebula 
or, in the case of the 1998 data, Mrk 421. Table | hsts the cuts used for each years' analyses. 



3.2. 7-ray Rate and Flux Calculation 

There are slightly different analysis methods for the on/off and tracking observation 
modes. For the on/off observations the background is estimated from the off source run, which 
is assumed to be on a sky region which does not include a 7-ray source. This analysis mode has 



been discussed at length elsewhere (e.g. Kerrick et al. 1995; Catanese et al. 1998). 



For TRACKING observations the background is estimated from the ON source run itself. All of 
the 7-ray selection criteria apart from orientation (a) are applied to the data. The background is 
then estimated from events which are not oriented towards the source. In this analysis, background 
events with values of a between 20° and 65° are used. Images having values of a between 65° and 
90° are discarded because of possible systematic effects due to truncation at the camera's edge. 
Once the number of events with orientations in the 20° to 65° range is known then the expected 
number of background events in the signal domain {a of 0° to 15°, or to 10° for data taken in 
1998) can be estimated. Off source data recorded for this source and others can be combined to 
calculate a ratio, r it Ar, of the number of events in the signal region to those in the 20° to 65° 
region in the absence of a source. In this case, the significance of a 7-ray excess, S*, is given by: 

/Non + rmoff + {Ar)m^f^ 
and the 7-ray rate {R ± AR) is calculated from: 



R±AR= — 2LL ± 1 -LL (2) 
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where Non is the number of counts in the 7-ray domain {a < 10° or 15°), Nofj is the number of 
counts in the 20° to 65° a range and t is the duration of the observation. 

The inclusion of the statistical error on the tracking ratio effectively limits the amount 
(duration) of tracking data which can be usefully analysed for an excess. Once the duration of 
the TRACKING data exceeds that of the off source data used in the calculation of the tracking 
ratio the error on the tracking ratio starts to dominate and to limit the significance of a detection. 

For the purpose of investigating possible differences in results produced by the on/off and 
TRACKING analyses, the Crab Nebula data were analysed with both methods ( iQuinn 19971) . The 
results demonstrated that the 7-ray rates derived using both methods were consistent and stable; 
the tracking analysis did not introduce any apparent variability in the rate. For the analysis 
presented in this paper, all of the ON source data were combined with the tracking data and the 
resulting database analysed using the tracking analysis. 

The data presented here were taken with different telescope configurations having different 
sensitivities. It is therefore necessary to normalize when comparing the different data sets. To 
a first approximation, this can be achieved by converting the rates to fractions of the 7-ray rate 
from the Crab Nebula taken with the same telescope configuration and analysed with the same 
cuts. The Crab Nebula is believed to be a steady source of VHE 7-rays, as has been observed 
by the Whipple Collaboration over the past decade ( Cawley et al. 1999| ). To convert a given 



7-ray rate to an integral flux, the rate as a fraction of the Crab Nebula flux was multiplied by 
(1.05 lb 0.24) X 10^^*^ cm~^s~^, which represents the integral Crab Nebula flux above 350 GeV 



(Hillas et al. 1998), the threshold of the analysis presented here. 



3.3. Test for Variability 

To search for temporal variability in the 7-ray flux we apply a test for a constant rate. 
The sum is converted into a probability (-Py-a) that the emission is constant about the mean 



using the incomplete gamma function gammq{a,,'x) ( Press et al. 1988 ) 



P^2=gammq\ — ^ — ,— \ (3) 

where N-1 is the number of degrees of freedom. The number of trials is taken into account by 
calculating the probability {Ptrials) of -f^^ occurring in N trials from: 

Ptrials = 1 - (1 - P^^f. (4) 



This method is used to test whether the distribution of measured 7-ray rates is consistent with 
statistical fluctuations about the mean for a range of time-scales. 



-7- 



4. Results 

The data have been analysed using the tracking analysis described above. This differs from 
that applied to the 1995 data by Quinn et al. (1996), in that a statistical error on the tracking 
ratio is now included and that a 10% systematic error is no longer added to TRACKING results as a 
careful study showed that the results of the ON / off and tracking analysis methods are in close 
agreement ( [Quinn 1997 ). In fact, these changes tend to cancel each other. 



For the 1995 data set as a whole we obtain a 7-ray rate of approximately 10% of that of the 
Crab Nebula (0.18 it 0.02 min^^, giving a 9.1cr excess). This rose to approximately 20% of the 
Crab Nebula rate for the following season (analysis of 1996 data reveals an excess of 11. la and a 
7-ray rate of 0.26 it 0.02 min~^), indicating that the average emission level had doubled since the 
previous year. 

The monthly and nightly average rates, in fractions of the Crab Nebula rate, over all four 
years of Mrk 501 observation are shown in Figure |. The rate in 1995 appears to have been 
constant with the exception of one night, MJD 49920, when the it was approximately 4. 60" above 
the average. A test gives a chance probability of 1.2 x 10"'^ (after accounting for trials) that 
the daily averages were constant during that month. The x'^ probability that the daily averages 
are constant over the entire 5 months of observation is 0.17 (after trials). The probability that the 
emission is constant when averaged on monthly time-scales is 0.06 (after trials). Conversely, in the 
1996 data set there were no obvious flaring episodes, but when averaged on time-scales of a month 
there is significant variability. The probability that the emission is constant for the monthly 
averages is 3.8 x 10~^, after accounting for trials. When each month is examined for variability 
with the rates averaged on time-scales of a day no significant variability is found. However, the 
X^ probabilities are smaller than for a similar analysis of the 1995 data, suggestive of increased 
day-scale flickering in the 7-ray emission. The results of the x^ test for variability of the monthly 
and daily averages are shown in Tables § and § respectively. 

Observations of Mrk 501 in 1997 began in January using the LZA technique ( [Krennrich et al. 



1998). The initial observations suggested that the flux was much higher than in previous years. 
Conventional SZA observations commenced as soon as possible in February. Other ground-based 
7-ray experiments were notified and high emission levels were verified. In March a joint lAU 
circular by the CAT, HEGRA and Whipple groups ( |Breslin et al. 1997| ) announced preliminary 
results. The Whipple observations continued through June and included observations made with 



the LZA technique. The results of the LZA data are presented elsewhere ( [Krennrich et al. 1998 ) 



Our SZA observations revealed that the VHE 7-ray emission was very strong throughout 
Spring - Summer 1997. Significant variability was observed in the monthly averages. The emission 
appeared to increase steadily from February through May and then to level off in June. The daily 
rates also exhibit dramatic variability. Day-to-day changes in the flux by factors >4 were observed 
and on eight occasions the flux more than doubled between consecutive nights. On four occasions 
there were equally rapid decays in the rate. The average flux level for the season is 1.4 times that 
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of the Crab Nebula - an increase by a factor of 14 from the level in 1995. The peak rate, observed 
on MJD 50554, is 3.7 times the Crab rate. The average rates for the 1997 data were calculated 
from only one run (the first with elevation above 55°) per night. This was done to remove a bias 
in the observing strategy whereby observations of Mrk 501 on a given night continued only if the 
source was in a very active state. 

The strong signal-to-background level in the 1997 data allowed a search to be made for 
variability on time-scales shorter than one day. For this test, data from each night on which there 
were three or more runs (approx. 1.5 hours, see §||) were analysed to test for run-to-run variability. 
A total of 24 nights satisfied this criterion. The probability that the emission was constant 
was calculated for each of the nights; the distribution of these probabilities is shown in Figure |2|a. 
For a statistically-variable source this distribution should be flat but there is an excess of nights 
with small probabilities. Assuming that the variations are purely statistical, the probability of 
getting seven nights in this first bin (width = 0.025) out of 24 trials is 1.45x10"^. Of these seven 
nights, two show statistically significant variations within themselves. The probability for constant 
emission on MJD 50577 is 5.2 x 10"^ while for MJD 50607 the probability is 5.8 x W'^ (after 
accounting for trials). The flux on MJD 50607 has a doubling time of ~2 hours. We thus identify 
these two nights as exhibiting significant hour-scale variability while the five other nights exhibit 
marginal variability. Figure ^ shows the 7-ray rates for the two nights with significant variability. 

A search of the 1997 data for variability on time-scales of less than half an hour has also 
been performed. For this test each of the 28 minute runs, 144 in total, was divided into three 
equal length intervals. Each triplet was then analysed for variability. No significant variations 
were found and the distribution of probabilities (Figure ^3) does not indicate any excess of low 
probabilities. Hence, we see no evidence for significant sub-hour scale variability. 

The average 7-ray rate for the 1998 data set was 0.42ib0.04 min^^, approximately 20% of the 
rate obtained from the Crab Nebula i.e. on average the emission was much lower than in 1997. 
In fact, the average rates for March, April and May are comparable to the level of the initial 
detection in 1995. There were however two significant flaring events. The first occurred in early 
March where an apparent rise and decay were observed. Unfortunately this flare is poorly sampled 
due to bad weather. The flux was observed to be relatively high (~1.3 times the Crab Nebula 
flux) on MJD 50876 and on the following night a flux of ~5.0 times that of the Crab Nebula was 
recorded. This is the largest flux detected to date from Mrk 501 by the Whipple Observatory 
10 m telescope. For the next observation on MJD 50880 the measured flux was still relatively 
high (~1.3 times the level of the Crab Nebula flux). A second flare occurred in June. The flux 
increased on two consecutive nights, peaking at ~1.1 times the Crab Nebula flux on MJD 50991, 
then decayed on a similar time-scale. The average flux level immediately to either side of this flare 
was below the sensitivity of the telescope. 
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5. Discussion 

We have demonstrated that rapid variabihty, a common characteristic of blazars at all 
observed energies, is also present in the VHE 7-ray emission from Mrk 501. Our 4 year data 
set spans a remarkable change in the flux level from Mrk 501: the average yearly emission level 
exhibited a fourteen-fold increase between 1995 and 1997 and the average daily flux varied by a 
factor of ~50 (see Figure ||). In 1997 large amplitude day-scale flares occurred frequently and were 
usually followed by equally rapid decays. Day-scale changes in flux by factors as large as 4.7 were 
observed. Episodes of significant hour-scale variability were detected, with one having a doubling 
time of 2 hours. In addition, there is evidence of consistent hour-scale variability which is not 
resolved in individual episodes. A variable flux from Mrk 501 was reported by at least four other 
atmospheric Cerenkov observatories in 1997 (see, e.g. Protheroe et al. 1998| ). 



The data presented here suggest that Mrk 501 was more variable when the flux level was 
higher. However, this effect could also be due to the sensitivity of the telescope. At low flux levels 
it takes longer to accumulate a significant signal, so the search for short-term variability may 
be limited by poor statistics. To address this issue we performed a test to see if the day-scale 
variability observed in 1997 would have been detected in 1996 and/or 1995 and if the month-scale 
variability seen in 1996 and 1997 would have been detected in 1995. We calculated the percentage 
deviations about the mean level from a period where significant variability was observed and 
then, using the mean signal and background level from another period, calculated the signal (and 
statistical error) which would have been observed given the same percentage deviations about that 
mean. The results are that neither the 1997 month-scale nor day-scale variability would have been 
detected at a significant level if present in 1995. However, the 1996 month-scale variability would 
have been significant if present in 1995 (chance probability of ~ 10^^), while the 1997 degree of 
day-scale variability, if present in 1996, would have been detectable (chance probability ~ 10"''). 
Thus we conclude that there was a change in the flaring characteristics, in addition to the change 
in the mean flux level between the different observing seasons. 

The VHE 7-ray emission from Mrk 501 exhibits rapid variability similar to that of the 
emission from Mrk 421. A major difference between the two objects is that the emission from 
Mrk 501 seems to have a base-level, which changes on monthly and yearly time-scales, whearas 
the VHE 7-ray emission from Mrk 421 has been described as consisting of a series of rapid flares 
with no underlying baseline ( [Buckley et al. 1996 ). The variability of the VHE 7-ray emission of 



Mrk 501 in 1995 and 1996 was similar to that at other wavelengths, with small amplitude, slow 
variations being more common than fast, large amplitude flares. 

The increase in the VHE 7-ray power and variability in 1997 was accompanied by an increase 
in the hard X-ray power and an extension of the synchrotron spectrum to at least 100 keV 



( Catanese et al. 1997 ). This is consistent with an inverse-Compton mechanism for producing 
VHE 7-rays. However, our observations of Mrk 501 to date cannot discriminate between the 
inverse-Compton (electron) models and those where the dominant particles producing 7-rays in 
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the jet axe protons. More densely sampled light curves, covering as broad a wave-band as possible, 
are needed to provide more insight into the mechanisms responsible for the VHE radiation from 
BL Lac objects. Future, more sensitive 7-ray observations of AGN with proposed detectors such 
as GLAST, HESS, MAGIC and VERITAS will allow the structure of flares on shorter time-scales 
to be determined. 
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Table 1. Summary of observations (combined hours of ON and TRACKING data) 





1995 


1996 


1997 


1998 


Feb 






3.3 


0.9 


Mar 


3.7 


18.0 


12.4 


5.5 


Apr 


17.7 


15.5 


20.7 


4.6 


May 


20.3 


10.3 


21.2 


20.5 


Jun 


11.5 


8.8 


9.7 


11.6 


Jul 


8.3 






8.7 



Total 61.5 52.6 67.3 51.8 
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Table 2. 7-ray selection critera 



Data set 


Uuts 




0.073° < width < 0.15° 




0.16° < length < 0.30° 


1995 


0.51° < distance < 1.0° 


and 


a < 15° 


1996 


size > 400 d.c. 




maxl > 100 d.c. 




max2 > 80 d.c. 




0.073° < width < 0.16° 




0.16° < length < 0.33° 




0.51° < distance < 1.17° 


1997 


a < 15° 




size > d.c. 




maxl > 95 d.c. 




max2 > 45 d.c. 




0.073° < width < 0.16° 




0.16° < length < 0.44° 




0.51° < distance < 1.25° 


1998 


Q < 10° 




size > d.c. 




maxl > 75 d.c. 




max2 > 65 d.c. 




nsym meln/ > 
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Table 3. Results of analysis for variability of the monthly averages 



Year Number of Months Rate ^ ^trials ^ 

1995 5 O.lOzbO.Ol 0.05 

1996 4 0.20±0.02 3.6x10-^ 

1997 5 1.39±0.07 S.SxlQ-^s 

1998 6 0.32±0.04 8.5x10-^^ 



'^Average rate over the observation period expressed 
as a fraction of the rate from the Crab Nebula. 

probability for constant emission after 
accounting for trials. 
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Table 4. Results of analysis for variability of the daily averages 



1995 



Obs. Period 


Number of Nights 


Rate ^ 


P , b 

^ trials 


Mar 28 


- Apr 07 


09 


0.09±0.02 


1.00 


Apr 22 ■ 


- May 10 


15 


0.08±0.02 


1.00 


May 21 


- Jun 05 


14 


0.07±0.02 


1.00 


Jun 18 - 


- Jul 01 


10 


0.10±0.03 


0.99 


Jul 22 - 


Jul 30 


08 


0.16±0.02 


1.8x10-3 


Total 




56 


O.lOibO.Ol 


0.22 


1996 


Mar 17 


- Mar 30 


13 


0.19±0.02 


0.33 


Apr 13 ■ 


- Apr 28 


15 


0.12±0.02 


1.00 


May 09 


- May 25 


15 


0.24±0.03 


0.18 


Jun 06 ■ 


■ Jun 21 


10 


0.29±0.03 


0.46 


Total 




53 


0.26±0.02 


3.8x10-6 


1997 


Feb 10 - 


■ Feb 16 


05 


0.62±0.08 


0.33 


Mar 05 


- Mar 16 


12 


0.96±0.08 


< 10-38 


Apr 07 ■ 


- Apr 19 


12 


1.37±0.08 


< 10-38 


Apr 30 ■ 


- May 15 


12 


1.68±0.09 


< 10-38 


May 27 


- Jun 11 


14 


1.61±0.09 


< 10-38 


Total 




55 


1.29±0.06 


< 10-38 


1998 


Feb 28 - 


• Mar 09 


06 


1.94±0.29 


< 10-38 


Mar 23 


- Mar 31 


03 


0.00±0.07 


1.00 


Apr 19 ■ 


- Apr 29 


09 


0.08±0.04 


1.00 


May 20 


- May 31 


10 


0.09±0.02 


0.13 


Jun 16 - 


- Jul 03 


15 


0.30±0.05 


< 10-38 


Jul 15 - 


Jul 27 


06 


0.12±0.04 


0.94 


Total 




49 


0.32±0.04 


< 10-38 



^Average rate over the observation period expressed as a 
fraction of the rate from the Crab Nebula. 

probability for constant emission after accounting for trials. 
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'a) Month ly Averages 
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Fig. 1. — Average 7-ray rates on (a) monthly and (b) daily time-scales for Mrk 501 between 1995 
and 1998. 
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- (a) Variability within eacin night 
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Fig. 2. — Distribution of probabilities for constant emission for (a) each of 24 nights (28 min. 
bins) and (b) each of 144 runs (9 min. bins) for data taken in 1997. There is a significant excess of 
small probabilities in plot (a), indicating hour-scale variability. 



-19- 



o 



o 
a: 



1 6 
1 4 
12 



6 
4 
2 




^ 14 
I 12 




8 
6 
4 
2 




'a) MJD 505: 



1 ^ 



9 



10 11 

UTC (hrs) 



(b) MJD 50607 




6 



7 



9 



10 11 

UTC (hrs) 



Fig. 3. — 7-ray rates for the two nights in 1997 which show significant variabihty on a time-scale 
of hours. The probabihty (after accounting for trials) for constant emission during May 9 (a) is 
5.2 X 10~^ while for June 8 (b) it is 5.8 x 10~^. The flux on the latter night has a doubling time 
of ^2 hours. 



